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Density of vortices and monopoles after a rapid 
pressure quench is predicted. Critical exponents are 
worked out in overdamped and underdamped limits. 



I. INTRODUCTION. 

Topological defects play a prominent role in many con- 
densed matter systems, see e.g. jl| for a review. It was 
also suggested that they were an important ingredient of 
the early universe S. Topological defects can be gener- 
ated in large numbers during a second order phase tran- 
sition. The dynamics of such a transition to a symmetry 
broken phase has been an object of much recent attention 
because of its importance in the cosmological context ||] 
and in condensed matter physics Q. An early estimate 
of defect density after a quench was given by Kibble §. 
It this theory, the speed of light is a dominant factor 
which determines the size of correlated domains. A more 
detailed scenario was put forward by Zurek j^] , who em- 
phasized the importance of the nonequilibrium dynamics 
of the order parameter. Recent experiments point to the 
latter theory. 

In a recent paper we calculated density of kinks in 
an underdamped one dimensional 4 theory. In this pa- 
per we give a comprehensive treatment of hedgehogs in 
D = 2, 3 spatial dimensions. A hedgehog is a point like 
defect which is topologically stable thanks to its asymp- 
tote at infinity. The asymptote is a topologically nontriv- 
ial map S D — > S D . A hedgehog should be distinguished 
from a texture the asymptote of which is a constant. The 
paper has essentially two parts. In the first we give a 
general outline of the formalism. In the second part the 
formalism is applied to vortices and monopoles. 



II. MODEL. 

Let us consider an O(D) symmetric time dependent 
Ginzburg-Landau model in D spatial dimensions 



(f) a (t,x) +j(f> a (t,x) 



d d 



where k = 1,...,D are spatial indices and a, = 
1, . . . , D are internal indices. Repeated indices imply 
summation. 4> a (t, x)'s are real. rj a (t,xys are Gaussian 
white noises of temperature T with correlations 

< ri a {t,x) >= , 
<V ai (h,x 1 )r ] a2 (t2,X2) > = 

2jT S aia2 5(h - t 2 ) 6^ {xt - x 2 ) . (2) 



The coefficient a(t) is time dependent, 
linear quench 
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< t < T^ 1 
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Before the quench, for t < 0, the system is in a symmetric 
phase (A > 0), during the quench, at t = t, it undergoes 
a transition from the symmetric phase a(t < r) > to a 
broken symmetry phase a(t > t) < 0. Finally it settles 
down at a(t) = —1. 



III. SOLITON TRANSFORM. 

The model ([!]) in its symmetry broken phase (a(t) 
-1) admits pointlike topological defects or solitons 



"(f) = S a {x-z). 



(4) 



z is an arbitrary location of the soliton. The soliton so- 
lution has D zero modes 



d 



S a (x-z) 







dx k 



: S a {x-z) . (5) 



cp a {t,x) 



dx k dx k 

-2a(t)<p a (t,x)-2[^(t,x)^(t,x)}<j, a (t,x)+7 1 a (t,x) , (1) <f> a {t,x) = S a [x~ z(t)} + 



When field equations are linearized around the soliton 
solution (Q), one obtains a discrete/continous spectrum 
of eigenmodes of the fluctuation operator. There are D 
zero modes ((^) and a spectrum of modes with positive 
eigenvalues and eigenfunctions P" (x — z) labelled by a 
discrete/continuous index p. A field <fi a (t,x) can be ex- 
panded as a sum of the soliton solution (^) with an appro- 
priate choice of the soliton position z and a combination 
of positive eigenvalue eigenmodes with coefficients $ p , 



j2* P (t)p?[z-m] ■ (6) 

p#0 
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p 7^ in the sum stresses that the zero modes are ex- 
cluded from the expansion, they are absorbed into the 
soliton location z{t) . Given the order parameter 4> a (t, x) , 
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one can trace the soliton position z(t) with a help of a 
" soliton transform" of the order parameter 



I k (t,x)= / d u yZ%(y)cl> a (t,x + y) 



(7) 



I k [t,z(t)] = for the field ^ because the zero modes 
(H) are orthogonal to the soliton solution (Q) and to the 
eigenmodes P" (x — z). 

The soliton position z(t) is a zero of the soliton tran- 

form, I k [t,z(t)\ = 0, §. 

It may happen that there are several solitons in the 
system. In that case an expansion like @ is possible 
around any soliton position Z\(t), &2(t), . . .. Each z(t) is 
a zero of the soliton transform. To count the solitons 
one should count the zeros of the soliton transform. An 
average density of solitons can be found as 



i(t) = 



d D x 
D L D 



<6^[I(t,x)} | 



d{x\... lX D ) 



(t, x) I > 
(8) 



where the integration goes over a D-dimentional box of 
size L. Following similar steps as in we obtain 

n(t) = Cd( 7- 1 r^) 2 , 

v ' v W{t,r = 0) ' ' 

W(t,\x\)=<I k (t,0)I k (t,x)> , 

if D = 1 

Co ■■■{ 4= ,if D = 2 (9) 
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Once the correlation function W(t, r) is known one can 
work out the average density of solitons in the system. 
In the considered cases the soliton transform does not 
distinguish solitons from antisolitons so Eq.(^) gives a 
total density of topological defects. 

It is convenient to express W(t, r) by a correlation 
function of the order parameter. The order parameter 
can be expanded in Fourier modes 



<P a (t,x)= / ^pe^f(t,p) 



(10) 



Thanks to translational invariance and the O(D) symme- 
try of the model (|l]), a correlation function of the Fourier 
modes has the form 



<r i (i,Pi)0 a2 (i,P^) > 



_ T6 a ^S<- D ^(p 1 -p 2 )G(t,p 1 ) 



(27T) 



D 



(11) 



This definition and the Fourier transform of the zero 
mode, 



(12) 



Ztiv)= J d D xe-^ s Z%(x) , 
U(\p\) = ZZ(+p)ZZ(-p) , 



leads to the following expression for the density of soli- 
tons 



i(t) = C D ( 



fd D PP 2 U( P ) G{t,p) o 



Jd D pU( P ) G(t,p) 



(13) 



To find the density of solitons just after the quench (||) 
we need a correlation function G(r,p) when the param- 
eter a(t) passes through the critical point and enters the 
symmetry broken phase. This correlation function is a 
topic of the next section. 



IV. CORRELATION FUNCTION AT THE 
CRITICAL POINT. 

As long as the system is still in the symmetric phase 
(t < r), the order parameter (f> a can be regarded as a 
small fluctuation around its symmetric ground state < 
<j> a (t, x) >= 0. It is justified to neglect the cubic term on 
the RHS of Eq. (|l|) . A Fourier transform of the linearized 
Eq.® is 

^\t,d + 4 a (t,pl = -lMt)+p 2 }4> a (t,p) + v a (t,p) , (14) 
<T(t,P) >=o , 
<v ai (h,Pi)v a2 (h,P2) >= 

2 7 T 



5 a ^ S(h - t 2 ) S^Hp! - p 2 ) 



(15) 



The correlation function G(r,p) can be easily worked 
out in two opposite limits. 



A. Underdamped limit, 7 



0. 



The underdamped limit is achieved if 2a(t) >> 7 2 /4 
from the beginning of the quench to the freeze-in time 
t. In this limit the system ( |l4| ) can be mapped to the 
equivalent model 



i (t,p) = -l4> a (t,p) + n a (t,p) , 
< n a {t,p) >= , 
<n Ql (ti,PiK 2 (t 2 ,p 2 ) > = 

2 7 T 5 aia2 Sjh-h) S^Hpx -P2) 



{2n) D [2a(t 1 )+pl} 



(16) 



To see the nature of this equivalence let us consider three 
examples. If a(t) = const, then both systems have the 
same correlation 



G(t,p) 



2a +p 2 



(17) 



Their response to the change of parameters is also the 
same. In particular they have the same relaxation time. 
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For example, if a(t) = const and the temperature jumps 
from T = for t < to T > for t > 0, then the 
correlation function for t > is 



G(t,p) 



1 



r[l-e-n 



2a + p 2 1 
2a + p^ 



-7tQ( 



\/2a+P 2 



)] (18) 



for the models (jig ) and ( [l4j ) respectively. 0(. . .) are os- 
cillating terms which are negligible for small 7. In this 
limit the response of the two systems is the same. An- 
other example is a response to a sudden change in the 
parameter a from a = a_ for t < to a = a + for t > 



G(t,p) 
G(f,p) 



1 



--ft 



-jt 



2a+ + p 2 
1 - er^ 



2a- 



-7* 



2a+ -J- 
7e 



0( 



p^ 2a_ 

-7T 



7e 



-7T 



^/2a+ + p 2 ' V 2 «- + J> 2 



(19) 



where 0(. . .) are once again oscillations negligible for 
7^0. The two considered systems respond in the same 
way to the impulsive perturbations of their parameters. 
The impulsive perturbations contain all frequencies so 
the models can be considered equivalent in the under- 
damped limit. 

In the underdampcd limit the relaxation times of all 
the Fourier modes are the same and equal to 1/7. One 
can talk about a relaxation time of the system as a whole. 
As a(t) begins to change at t = Q the system tries to 
adapt to the changes. It can not react faster than on the 
time scale of its relaxation time 1/7. The correlations 
are effectively frozen at the freeze-in time i, when the 
time left to the transition is equal to the relaxation time, 
(t — t) = -. At this instant a(t) — so the frozen 
correlation function is 



1 



M + p 2 



The frozen correlation length can be identified as 



-1/2 



(20) 



(21) 



These asymptotic results are valid for r >> 1/7. The 
system remains in the underdamped limit until the 
freeze-in time, 2a(t < t) >> 7 2 /4, if the condition 



8A > T7 3 



(22) 



is satisfied. If r << I/7, then the system is frozen from 
the beginning of the quench and the frozen correlation 
length is |= 1/V2A. 

The analysis so far was on a heuristic level. The equa- 
tion (|l6|) is of first order in time, so the exact correlation 



function of the system ([L6|) can be easily worked out. The 
formal solution of the model ( |l6| ) is 

4> a (r,P) = C dti e^-WrWtup). (23) 



(24) 



With this solution the correlation function is 



G(r,p) = 



~7T 



2A + p 2 



dx 



^x+p 2 

7T " 



It turns out that Eq. (|2(]) is indeed the asymptote of the 
exact solution for r >> 1/7. 



B. Overdamped limit, 7 — > 00. 

In the overdamped limit (7 2 /4 >> 2a(t)) the second 
order time derivative on the LHS of Eqs.(|jl4) can be 
neglected. Eq.(jlJ) is simplified to 



7^ (t, p) = - [2a(t) + p 2 W (t, p) + r (t, p) 



(25) 



In an overdamped system a relaxation time of a given 
Fourier mode depends on its momentum p and equals 
7/2[2a(i) +p 2 ]- There is no universal relaxation time of 
the system as a whole. Short wave lengths freeze in later 
than the long wave lengths. For a given momentum p, the 
relaxation time is equal to the time left to the transition, 
if (t — t p ) = j/2[2a(t p ) + p 2 ). On this heuristic level the 
long wave length limit of the frozen correlation function 
is 



G 



heuristic 



( T >p) 



1 



2A +p 2 



The frozen correlation length can be identified as 



-1/4 



^heuristic 

The model (^5|) has an exact solution 
~^{T,p) -- 



(26) 



(27) 



dh exp{- / dt 2 [2a(t 2 )+p 2 }/^}v{t u p) 

(28) 



This solution, the correlations in Eq. ([14|) and the explicit 
form of a(t), see Eq.(||), give an exact correlation function 
at t = r 



G(t, P ) = 



e -2p 2 r/ 7 

2A+p 2 



2-fA 



2-fA 



(2A + p 2 ))-Erf{ 



2-fA 



(P 2 ))} ■ (29) 



If t is much longer than the relaxation time of the p = 
mode before the quench, r >> 27^4, the correlation 
function takes an asymptotic form 
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G(t, P ) « J— e*^ / dxe" 21 . (30) 

The frozen correlation length can be identified as 

i^o = (^)* • (3D 



V. DENSITY OF VORTICES. 

In two dimensions (D — 2) the model (0) admits a 
topologically stable vortex solution 

S a (x) = f(p){6 al cos(<p)+5 a2 sm(<p)} , (32) 

where p, ip are polar coordinates and the profile function 
interpolates between /(0) = and /(oo) = 1. A single 
vortex has infinite energy, the IR divergence is due to the 
long distance asymptote 



S«^) = r 1 cos(^)+^ 2 sm^)] 



(33) 



The asymptote is scale invariant. The asymptote de- 
termines the long wave length asymptote of the Fourier 
transform, which is 



U(p) 



up to a numerical factor. 



1 



(34) 



A. Underdamped limit. 

Combination of the formulas (|2(||l3],[34|) should give the 
density of vortices for large r. The numerator integral in 
(O) is proportional to 



+ OC 



p dpp 2 U(p) G(r,p) 



-f-oo 



p dp p U (p) 



2A +p 2 



p dp U(p) A w° - ln(A) 



(35) 



A is an IR regulator. The denominator integral is pro- 
portional to 



-f-oo 



p dp U(p) G(r,p) 



A 



+ °° A TTt ^ 1 



2AJ A p 



2A 



We manipulated with the order of the limits r — > oo and 
A — > so as to get rid of the IR divergence and obtain a 
finite regularized asymptote 



n REG. 



A 

7T7T 



(37) 



The results for vortices are dominated by the IR diver- 
gence, which is due to the asymptote (|33|). The asymp- 
tote is scale invariant so the density of vortices is just an 
inverse squared of the frozen correlation length (pi]) . 



B. Overdamped limit. 

In this limit we combine the formulas ( |30| , p^| , p^ ) . The 
numerator in Eq.(|l3|) is proportional to 



p dp p 2 U (p) G(r, p) 



+ OC 



pdpU{p) A «°-ln(A) 



(38) 



The denominator is 

TO 

p dpU(p) G(r,p) 



-Hoc 



/ 2r / j th \ ze- I j -x 2 

' — - / pdpU{p)e 2 ~< A I dx 6 ~ 



(39) 



The regularized asymptote of the density of vortices is 

1 



n REG. 



rfA 

27r 3 / 2 V ~ 



(40) 



The 1/2 critical exponent is consistent with the numerical 
results in pCfl . 



VI. DENSITY OF MONOPOLES. 

The analysis of monopoles in D = 3 dimensions is 
much the same as that of vortices in D — 2 dimensions. 
The monopole configuration is 

S a (x) = f(r)[S al sin(0) cos(<^) + S a2 sm(6) sm(ip) + S a3 cos(0)] 

(41) 

where r, 9, ip are spherical coordinates. The profile func- 
tion interpolates between /(0) = and /(oo) = 1. A 
single monopole has infinite energy, the IR divergence is 
due to the long distance asymptote 



(36) % {x) = [S ai sin(0) cos(y>) + S a2 sin(0) sin(^) + 8 a6 cos(0)] 

(42) 
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The asymptote is scale invariant. The asymptote deter- 
mines the long distance behaviour of the Fourier trans- 
form, which is 



U(p) 



1 



(43) 



up to a numerical factor. Similarly as for vortices this 
singularity alone determines the critical exponents. 



A. Underdamped limit. 

Combination of the formulas ( pfj| , [T^ , ^3| ) gives the den- 
sity of monopoles for large r. The numerator integral in 
( |l3| ) is proportional to 



p 2 dpp 2 U(p) G(t, P ) 

' 1 
p 2 dp p 2 U (p) 



2A + 2 
■yr r 



" a a^o 1 



p dp U(p) 



A 



(44) 



A is an IR regulator which is provided either by the sys- 
tem size or by an average distance between monopoles 
and antimonopoles. The denominator integral is propor- 
tional to 



+ OC 



p 2 dp U{p) G(t, p ) 



+ 30 



p 2 dp U (p) 



2A 
7T 



7T f dp 7T 1 
2AJ A p^~ 2AA 



(45) 



We manipulated with the order of the limits r —> oo and 
A — > so as to get rid of the IR divergence and obtain a 
finite regularized asymptote 



A(-> ! 



(46) 



The asymptote (£3|) is scale invariant; the density of 
monopoles scales like an inverse qubed of the frozen cor- 
relation length (pl|). 



B. Overdamped limit. 

In this limit we combine the formulas (|3(],[l3||l3|) . The 
numerator integral in Eq.(|l3|) is proportional to 

p 2 d PP 2 U(p) G(r,p) T ~°° 

A 

P d P U(p) ~ - . 

A A 



(47) 



The denominator is 




p 2 d P U(p)G(T,p) ; 
' J p 2 dp U(p) 



p dp U(p) 




dx e 



„2 A^O 



(48) 



The regularized r >> 27^4 asymptote of the density of 
vortices is 



nREG. 



—0—v 



(49) 



VII. CONCLUSION. 

We considered density of hedgehogs just after a rapid 
pressure quench. For large enough quench time r the 
density depends on r like 



n ~ r " . 
Our findings are summarized in a table. 



(50) 



D 


DEFECT 


a 


7^0 


7 — * oo 


2 


VORTEX 


1 


i 

—I— 


3 


MONOPOLE 


3 

2 


4 



The critical exponent depends on the degree of damp- 
ing. The underdamped critical exponent is twice the 
overdamped critical exponent for a given defect. It is not 
surprising as the two limits are dominated by a second 
order and a first order time derivative respectively. 

Global vortices and global monopoles interact by long 
distance confining potentials. Their density just after 
a quench is merely an initial condition for their subse- 
quent rapid recombination described by phase ordering 
kinetics [ HJ . Their critical exponents apply to local vor- 
tices/monopoles coupled to a gauge field provided that 
the gauge field equilibrates with the matter field rather 
slowly |1(| . These local defects have finite energy and in- 
teract weakly so their subsequent recombination is purely 
diffusive. 

Topological textures like Bloch waves for D = 1, 
skyrmions for D = 2, magnetic bubbles and Hopfions 
for D = 3 require a separate study |12|. 
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